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Abstract. This paper presents a study of the free energy and particle density of 
the relativistic Landau problem, and their relevance to the quantum Hall effect. 
We study first the zero temperature Casimir energy and fermion number for Dirac 
fields in a 2+1-dimensional Minkowski space-time, in the presence of a uniform 
magnetic field perpendicular to the spatial manifold. Then, we go to the finite- 
temperature problem, with a chemical potential, introduced as a uniform zero 
component of the gauge potential. By performing a Lorentz boost, we obtain 
Hall's conductivity in the case of crossed electric and magnetic fields. 
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1. Introduction 

The quantization of the Hall conductivity Tl is a remarkable quantum phenomenon, 
which occurs in two-dimensional electron systems, at low temperatures and strong 
perpendicular magnetic fields. Most proposed explanations for this phenomenon (for 
a review see, for instance, |2|) rely on one-particle theory and make use of the Kubo 
formula concerning the conductivity as a linear response function to the external field 

El . . . 

It is the aim of this paper to show that, in the context of relativistic field theory, 

the quantization of the Hall conductivity in multiples of the magnetic flux arises as a 
consequence of the spin-statistics theorem, which is a straightforward outcome of such 
a theory. 

After this calculation was finished, a very interesting, recently published |^, 
calculation of the Hall conductivity in relativistic systems was brought to our 
knowledge. We will compare our results to those obtained in this recent publication 
in our final comments. 

In section 121 we present the theory of Dirac fields in 2 + 1 Minkowski space-time, 
interacting with a magnetic background field perpendicular to the spatial plane, and 
evaluate the vacuum expectation values of the energy and fermion density. 
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Section |51 contains the generalities of the same theory in EucUdean 3-dimensional 
space, and presents the eigenvalues of the corresponding Dirac operator. From such 
eigenvalues, the partition function is evaluated in section^ 

Section |S1 contains the resulting free energy and mean particle density at finite 
temperature, and a study of their zero-temperature limit. In the same section, we 
perform an adequate Lorentz boost in order to consider the problem of fermions 
interacting with crossed electric and magnetic fields, and obtain the Hall conductivity. 

Some final comments about the case in which the chemical potential /x coincides 
with an energy level are presented in sectional In the same section, we discuss the 
role that experimental results on the quantum Hall effect can play, as an arena for 
testing the physical relevance of the phase of the determinant or, equivalently, of the 
multiplicative anomaly, in view of our discrepancies with reference ^ . 

A very sketchy presentation of the present study can be found in 0. 



2. Zero-temperature problem 

We will use the metric (— , +, +) , natural units Ti = c = 1, and choose the following 
representation for the Dirac matrices: namely, 

7a/ = *c^3, 7m = and 7^,^ = cti. (1) 

The Hamiltonian can be determined from the solutions of the Dirac equation 
eyl)^' = 0, where — e is the negative charge of the electron. In the Landau 
gauge A — (0,0, -Bx) , with B > Q. Thus, after setting ^>{t,x,y) = e~^^*'ip{x,y), we 
get the Hamiltonian H = icridx — i<J2dy + <J2eBx. The corresponding eigenvalues can 
be determined from 



—E idx — dy — ieBx 

idx + dy + ieBx ~E 



^{x,y) =0. (2) 



In order to solve the above equation, we take 



(Pk{x,y) A 1 / e^'^y^kix) 



This leads to the following system of first order equations 
{idx - ik - ieBx)xk = Etpk 

{idx + ik + ieBx)tpk = Exk ■ (4) 

We first solve for the zero mode Eq — 0. In this case, both equations decouple 
and, after imposing that the eigenfunctions be well-behaved for x — > ±00 we obtain, 
up to normalization 

_ f (eS/7r)i/4exp{~ieB(.T + fc/eS)2} 



Mx) = ^ ' ' 2--^-^-/--^ / j . (5) 

In the case E ^ 0, the normalized eigenfunctions which arc well-behaved for all 
values of x can be written in terms of the Hcrmite's polynomials Hn{u) 6 as follows 

^""^"^ = exp{-iei3(x + VeB)^} (6) 

±Hn {VeB{x + k/eB) 



iV2n H, 



n-i {x + k/eB)^ 
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with corresponding eigenvalues 



En — iV2neB , ?i — l,...,oo. (7) 

We notice that each eigenvalue exhibits the well known Landau's degeneracy per unit 
area: namely, 

zvr 

Had we chosen the other non-equivalent representation of the gamma matrices in 
2 + 1 dimensions, exactly the same spectrum would have been obtained obtained, the 
only difference being the chirality of eigenfunctions. 

The vacuum expectation value of the energy per unit area, defined through a zeta 
function regularization (see, for example, T*, and references therein), is given by 



(9) 



In the present case, we have (a is an arbitrary parameter with mass dimension, 
introduced to render the complex powers dimensionless) 



(10) 



s=-l 



Always in the zeta-function regularization framework, the fermion number is |S] 

\B„>0 E„<0 /J 

where Nq is the contribution coming from zero modes. 

In our case, the nonvanishing spectrum is symmetric. So, only the zero mode, 
which is charge self-conjugate, contributes. This gives as a result [S] 

N{B)=±^. (11) 
Or, equivalently, for the vacuum expectation value of the charge density 

j°(i?) = Te^. (12) 

The sign indetermination is a natural consequence of the twofold vacuum 
degeneracy. 



3. The theory at finite temperature with chemical potential 

In order to study the effect of temperature, we go to Euclidean space, with the metric 
(-I-, +, +) . To this end, we take the Euclidean gamma matrices to be 70 = = — ^3, 
7i = ^l^j = (72, 72 = 7^ = Ui- We will follow 9 in introducing the chemical potential 
as an imaginary ^0 = ^^■f i^i Euclidean space. Thus, the partition function in the 
grand-canonical ensemble is given by 

\nZ det{ip~ e^) . (13) 
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In order to evaluate the partition function in the zeta regularization approach 
| 1(J) . we first determine the eigenfunctions, and the corresponding eigenvalues, of the 
Dirac operator, in the same gauge used in the previous section, i.e, we solve 

[-ia^idr + ^J)+ i(J2dx + (Ji{idy - eBx)]^ = uj-^ , (14) 



or, after writing '^[t,x, y) 



-i{dr + ^J.) dx+idy — eBx\f $(t, a;,y) 



-dx + idy — eBx 



i{dr + fJ-) 



^{T,x,y) 



^ $(t,x,2/) 
2(r,a;,y) 



(15) 



In order to satisfy antiperiodic boundary conditions in the r direction, we propose 



'^kA{T,x,y) = 



with 



Xi = {2l + 1)-, 



(16) 
(17) 



where /3 — is the inverse temperature. 
After doing so, and writing 



Xk,iix) 



we have, for each k,l, 

{dx - k ~ eBx)xk.i =[t^-h)'Pk,i 

{-dx - k- eBx)(pk,i = (t^ + Xi)Xk,i , 

where A; = A; — ifi. 

There are two types of eigenvalues and corresponding eigenfunctions 



(18) 



1) uJi ~ A/, with / = — oo, ...,cxD, and corresponding normalized eigenfunctions 



= (f) 







(19) 



Note that these eigenvalues are not square roots of the eigenvalues of the squared 
operator. They will eventually lead to a "spectral asymmetry" || and, thus, to a phase 
of the determinant, which will be studied in detail in the next section. 



2) W!,Ti = ±Y'A^ + 2nei?, with n = l,...,cx), I = — oo,...,cxd, and corresponding 
normalized eigenfunctions 



i'k,l,n{x) = Ak,l,, 

where 



A _ , eB 

Ak,Ln — I 



2 — 



Try [{n-iy.]i 



2neB 



2neB+ |w/.„ + A/|2 



(20) 



II Here, quotation marks are due to the fact that, the Dirac operator not being self-adjoint, we have 
a complex spectrum. 
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In all cases, the degeneracy per unit area is again given by in equation 
Choosing the other nonequivalent representation of the gamma matrices leads to a 
change in the eigenvalues of type 1), which are replaced by uji = —A;. However, as 
will be discussed below, this doesn't lead to a change in our physical predictions. 

4. Evaluation of the partition function at finite temperature and chemical 
potential 

The partition function, in the zeta regularization scheme JOI j is given by 

d 



log 2 = 



ds 



C{s, — : — ) 



(21) 



As in the previous section, a is a parameter with mass dimension, introduced to render 
the (^-function dimensionless. 

We must consider two contributions to log Z, respectively coming from eigenvalues 
of type 1) and 2) in the previous section, i.e., 



where 



Ai(/i) = - 



ds 



s=0 



(2^ + 1)^-^ 

afJ a 



(22) 



(23) 



and 



where 



A2(/i,i3) 



d 
ds 



C,2{s,^i,B) , 



s=0 



C2(s,Ai,B) = (l + (-l)~^)AL 

n = l 
I = —oo 



+ l{2l + l)—-i^^ 
\ ap a 



(24) 



(25) 



The contribution Ai{fi) can be evaluated at once for the whole /i-range. The 
analytic extension of Ci(s,m) can be achieved as follows (for a similar calculation, see 

El) 



Ci(s,Ai) = Ai J2 



1 — — OQ 



27r 



a(3 

Ar. ( — 

ap 



E 

.1=0 



I — 

a 






*27r_ 


2n 



— s oo 



E 

(=0 



1=0 



— (/ H — ) — I — 
^ 2^ 2tt 



(26) 



Now, in order to write the second term as a Hurwitz zeta, we must relate the 
eigenvalues with negative real part to those with positive one without, in so doing. 
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going through zeros in the argument of the power. Otherwise stated, we must select 
a cut in the complex w plane [r2] . This requirement determines a definite value of 
(— i.e., (—1)^'* = ^nrsign(ti)s ^ Taking this into account, we finally have 



Ci(s,m) 



2lT 

13a 



+ e 



2tt ) 



(27) 



From this last expression, the contribution Ai(/i) to log 2 can be obtained. It is 
given by 



Ai(a*) 



0, 



1 



A, 



Al <^ log 2 cosh 



1 



ifif3 
'2^ 



2 



(28) 



As commented in advance, the other nonequivalent representation of the gamma 
matrices leads to the same result for this contribution. In fact, even though this part 
of the spectrum changes sign, such change is compensated by the selection of the cut 
in the w-plane. In this case, one has 



Ci(s,Ai) = Al (I^ 
which also leads to H28|) . 



^"^'2 2^ 



The analytic extension of C2(sj B) requires a separate consideration of different 
fi ranges. We study in detail three of these cases. The generalization to arbitrary 
/z-ranges will be evident from these results. 

4.1. ^2 < 2eB 



C2(s,m,S) = (1 + (-1)~^)Al 

n = 1 '- 

I — — oo 



2neB 



+ (21 + 1)— -i^ 



a/3 



a 



Making use of the Mellin transform, this can be written as 



C2(S,^,S) = 



(1 + (-!)-«) A, 



dtt^' 



oo 



- + ((2'+l)^-''^)' 



n = 1 

I — — oo 



or, equivalently, 



(29) 



(30) 



(1 + (-1)-)Al ^ 
C2(s, Ai, B) = 2^ 



X 63 



n=l 

~'2t / IT i/i 

a(3 \aP a 



dtt-^ 'e 

AlTt 



where we have used the definition of the Jacobi theta function 



63(2, :r)= e-^-'' 



(31) 



(32) 



/ — — 00 
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To proceed, we will use the inversion formula for the Jacobi function 

y/x \IX X J 

thus getting 

(1 + (-1)-'^) °° '■^ 

71=1 • 



C2(s,Ai,B) 



2 27r 47rt 
Applying once more the definition (|32|l . we have 

(f + (-f)-^)ALa/3 



2y^r(f) 



{/.oo OP 
/ dtt^^^j:, 



2/ dtt^-i ^ (-f)'cosh(/i/3Z)e 



After performing the integrals fQ, we obtain 



2v^r(f) 



1\ /2eB 



s - 1 



'SneB 



or, making the simple pole of the F function explicit 



C2(s,M,-S) 



(l + (-l)-^)ALa/?g 

4v^r(^) 



s - 1\ /2eS 



s - 1 



4 D-i) 



n,l=l 



cosh(^/3Z)i4:^ V2neB/32/2 



From this expression, the contribution A2 to the partition function can 
obtained, since the factor accompanying s is finite at s = 
-Al/3 ' 



A2(/i,S) 



20F 

00 



F(-i)\/2^C/?f-^ 



2o2 \ -4 



cosh(/x/3Z)i^„i W2neB(3H'^ 



n,l=l 

After using that 

1 (x) 

it can be written as 

A2(m,B) = Ai/3 



V2ZBCfi ( -- 



2 gLl)!cosh(M/?Oe-^^' 
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The sum over / can be explicitly performed, to obtain 

1' 



A2(Ai,S) = Ai/3 



R 



. oo 

i log (l + e-2V2^I5B/3 ^ 2 cosh (A*/3)e-^^^) 

n— 1 



(40) 



Finally, adding the contributions given by equations l|28f) and 1)41) |l we obtain, for 
the partition function in the range fj,^ < 2eB 



log Z = Al <^ log 2 cosh 



fi(3 



2 , f3V2^CR , 2 



+ 5^1og(l + e-2^^/3+2cosh(/x/3)e-^^'^ 



(41) 



4.2. 2eB < fi^ < AeB 
As before, we have 



2neB 



(2^ + 1)^-^ 

ap a 



C2(s,M,s) = (i + (-i)-^)Ai 

I = — oo 

However, in this range of /i, the contribution due to n = 1 is given by 

d 



ds 



cr\s,^^,B), 



s=0 



where 



(42) 



(43) 



cr'(s,A^,s) = (i + (-i)-^)Ai 



l — ~-oo _ 



2eB / . , . TT // 
-^+((2/ + l)^-z^ 



a/? a 



(44) 



The analytic extension of this expression must be performed in a different way. In 
fact, the expression cannot be written in terms of a unique Mellin transform, since its 
real part is not always positive (note, in connection with this that, for n — 1, equation 
H37|l diverges). Instead, it can be evaluated as follows 



1=0 



2eB+ ( {2l + ^)--i^l 



W2eB+(2l + l)--iiJ. 

P 



-iV2eB + {21 + 1) i/i 

P 



This can be written as a product of two Mellin transforms 



(45) 
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(i + (-ir^) 

a-[r(|)]2 



Jo 

or, after a change of variables 



(46) 



Now, calling x = |-, one has 



i=0 







-z^-t^' 


/ 


1 dt 








4 



2«-ir(§)]2 



2«-ir(§)]2 



^-[(2/+l)^-ip]2 



(47) 



oo 


■z2" 

4 _ 









(48) 



Finally, the x-integral and the sum of the geometric series can be performed to 
obtain 



-1 



dzz^Js-i{V2eBz)- 



+ (49) 



'sinh(f ) 

Now, the integral in this expression diverges at z = 0. In order to isolate such 
divergence, we add and subtract the first term in the series expansion of the Bessel 
function, thus getting the following two pieces 



X / rfz Z 2 

'o 



sinh(f ) 



and 



2^a-r(|)r(^)--yo sinh(f) 



dzz"-^ 



(50) 
(51) 



The contribution of equation (|50|) to the partition function, defined as in equation 
H43|l can be easily evaluated by noticing that the factor multiplying s is finite at s = 0. 
Thus, one has 



A^-7A(^,B) = -Al / dzz-^ cos(V2^^z)-l 
Jo L 



sinh(2^) 



(52) 
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where we have used that J_i{\/2eBz) - 
with fj, — !■ — /i, one can change z — > — z to obtain 

poo 

A^Y/)(^,B) = -Al / dzz~^ cos (\/2lBz) - 1 



,^ cos(V2ei3z). Now, in the term 



sinh(^) 



(53) 



This last integral is easy to evaluate in the complex plane, by carefully taking into 
account the sign of /i, as well as well as the fact that 2eB < /i^ in closing the integration 
path, to obtain 



1=1 



i-iy 



cosh(V2iB/30e 



(-1) 



i+i 



(54) 



MP 



or, after summing the series 

A^jl^{fi,B) = Al {log (l + e-^l'^l'^ + 2 cosh (72^5/3)6- 

+ |M|/?-21og(^2cosh(^)^| . (55) 

In order to obtain the contribution coming from H51|l . the integral can be evaluated 
for dis > 1, which gives 

(1 + (-l)-^)r(s)0F(a/3)^Ai 



(27r)"2«-ir(f)r(^) 

1 , ii-10 , 1 , iud 
Ch(s,x(1-— )) + Ch(s,;.(1 + — )) 



(56) 



where C,h(s, x) is the Hurwitz zeta function. The contribution to the partition function 
can now be evaluated by using that Cff (0, ^(1 - ^)) + Cj?(0, ^(1 + ^) = and the 
well known value of —-j^\s=qC,h{six) "5", to obtain 



A^,72')(/AS) = 2Ailog(2cosh(^)), 



(57) 



Summing up the contributions in equations H28|l . H55(l and H57|l . as well as the 
contribution coming from n > 2, evaluated as in the previous subsection, one gets for 
the partition function 



logZ = Al <^ log 2 cosh 



log(l + e~2|''l'' + 2cosh(\/2iB/3)e"l^l'^) + pV2^ ('''" ("^) " 

OO ^ 

^ log (l + e-2V2;^/3 + 2 cosh {^l(5)e 



-V2neBP 



(58) 



At first sight, this result looks different from the one corresponding to /i^ < 2eB 
(equation H41|l '). However, it is easy to see that both expressions coincide. The only 
difference is that (I58|) explicitly isolates the zero-temperature behavior from finite- 
temperature corrections for this range of fi. 
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4.3. ^2 ^ 2eB 

In this case, the analytical extension can be performed, exactly as in the previous 
subsection, up to equation (|52|l . Then, the term with /i —>■ — /i can be explicitly 
summed, which gives 

POO 

A^yAd^l = V2^,B) = -Al / dzz-'^ 

Jo 

This integral can be found in [H], and it gives as a result 

A^Jl^i\^^\^V2^,B)^AL 

The contribution A^^j^^{\n\ = V2aB, B) can be evaluated exactly as in the 
previous subsection. It is given by 



cos 



I cos Aeaz 
sinh {T^^ 



(59) 



log (^cosh [^1eB(i^ ^- 2 log | cosh 



\/2i^/3 



(60) 



A^J^VIa^I = V2^,S) = 2ALlog(2cosh(^-^)) . (61) 

These two contributions, together with those coming from n > 2 and from 
equation H28() . finally give for the partition function at this particular value of /i the 
same result as H41|l or H58II evaluated at ^ — 2eB. 

Thus, for any range of /i, one has 

logZ = A. {log (2 cosh (f )) - ^ + pV^Cn (4) 

00 ^ 

+ J2 log (1 + + 2 cosh (^/3)e-^^'^) , (62) 

which is continuous, even when /i coincides with an energy level. 

When this calculation was finished, it was pointed to us that, in reference the 
partition function for Dirac fields interacting not only with the background magnetic 
field, but also among themselves through a Nambu-Jona-Lasinio term, had been 
obtained through a different regularization (see also jj^; for a calculation of parity- 
breaking corrections, see JSl)- However, in the absence of this last interaction, the 
result in such reference differs from the present one for two reasons: In the first place, 
the authors of \1'6\ were considering a reducible representation of the gamma matrices, 
which lead them to an overall factor of 2 with respect to (|62|l . But, more important, in 
the same reference the determinant of the Dirac operator was evaluated as the square 
root of minus the squared Dirac operator, thus not including the factor coming from 
the phase of the determinant or, equivalently, from the multiplicative anomaly |16j . 
which we did consider (see equation (|2SJ)). As we will discuss in detail in sectional 
this leads to completely different predictions regarding the Hall conductivity. 
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5. Free energy and particle density 

From equation H62|l. the free energy per unit area [F = ^ogZ) can be obtained ^. 
It is given by 



F{y.,B,(3) - |i log (^2 cosh (^^)) - M + 

+ ^ E (l + e-^^^P + 2 cosh {txP)e-^^^) \ , (63) 

Moreover, the free energy is continuous at = 2keB,k — 0, ...,oo. In the low- 
temperature limit one has 



F{2keB < ii' <2{k + l)eB) -^p^^ -Al\^V2^ (^r.(^-^ " ^ j • 

Note that, for fc = 0, this result coincides with the Casimir energy obtained in 
section|21 even for /i 0, but in this range, i.e., for /j, less than the first Landau level, 
if positive, or greater than minus the first Landau level, if negative. 

The mean particle density can be obtained, also from (|62|l . as iV = i^logZ. 
For nonzero temperature and arbitrary /i one has 



N{^l,B,p)^^L\]- 



tanh (^^) ~ sign{^) 



^ 2sinh {fi(3)e-'^^^'^ 

2-^^ 1 + f,-2V2^p + 2 cosh {p(3)e^^^^^ 



(64) 



For nonvanishing /i, the low-temperature limit differs, depending on the /i-range 
considered 

N{2ekB < ^? < 2e{k + l)B) -^p^^ k/^L signin) , 
where, as before, k 



2eB 

This result was to be expected for particles with the statistic of ferniions, 
since relativistic field theory naturally leads to the spin-statistics theorem. At zero 
temperature, fj, is nothing but the Fermi energy; for example, for /i > 0, as ^ grows 
past a Landau level, such level becomes entirely filled. 

From the previous result, the mean value of the particle density at zero 
temperature can be obtained. After recovering units, one has 

—kce^B 

f{2ec^hBk < fi^ < 2eBc^h{k + 1)) = sign{fj.) , 

the other two components of the current density tri-vector being equal to zero in 
the absence of an electric field. Now, the zero-temperature limit of the same tri- 
vector in the presence of crossed homogeneous electric {F') and magnetic {B') fields 
can retrieved, for F' < cB' , by performing a Lorentz boost with absolute value of 
the velocity -gy. Suppose, for definiteness, that the homogeneous electric field points 
along the positive y axis. Then, the velocity of the Lorentz boost must point along 
the negative x-axis, and the transformation gives as a result 

.,0 -nce^B' . , . .,x -ne^F' . , . .,y „ 
J = sign{fi) , J ^ sign{^i) , j = . 

% Consistently with the comments in previous sections, all the results in this section and in the rest 
of this paper are independent from the representation of the gamma matrices chosen. 
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As a consequence, the quantized zero-temperature Hall conductivity is 

9 

— ne 

(Txy = sign{fi) . 

6. Final comments 

A comment is in order concerning the zero temperature value of the particle density 
in equation H64() when the chemical potential coincides with an energy level, i.e., for 
/Lt^ = 2ekB, k = 0, ...,00 In all these cases, the operation of setting /i to its value 
doesn't commute with that of taking the (3 —^ 00 limit. 

For instance, in the case n — 0, N{0, B, (3) is undefined at all temperatures, while 
lim^^o [lim^^oo iV(/i2 < 2eB,B,/3)] = 0. 

On the other hand, each time /i^ = 2ekB, k ~ l,...,oo, one has 
\iuYp^ooN{±V2ekB,B,P) = ±Al(A; - i), while \iin^^^^^[Yinii3^^N{ii,B,P)] 
is undetermined. 



As already commented in the Introduction, after this calculation was finished, 
reference ||4^ was brought to our attention. Their prediction concerning the Hall 
conductivity differ from ours. In fact, the Hall conductivity in that paper contains 
an overall factor of 4 with respect to ours, which is due to the use of a reducible 
representation and an extra sum over two values of spin. This is a more or less 
trivial difference. But, more important, after this factor is removed, their Hall 
conductivity is quantized in half-integer units of magnetic flux. This originates the 
word "unconventional" in the title of Our calculation, instead, leads to a Hall 
conductivity quantized in terms of integer quanta of flux density, which is entirely 
conventional, but obtained here from very first principles. The difference between 
both results is due to the inclusion, in our calculation, of the phase of the determinant 
(or, equivalently in this case, of the multiplicative anomaly). The interesting comment 
at this point is that experimental results on the Integer Quantum Hall Effect can 
clarify the physical relevance of the multiplicative anomaly (some recent results (T7II18) 
seemingly favor "unconventional" quantization). 

Finally, we mention that the more realistic case of massive fermions is at present 
under study [T!?] . 
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